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We obtain the electric field and scalar field for a static point charge in closed form in the 5D 
Schwarzschild-Tangherlini black hole spacetime. We then compute the static self-force in each of 
these cases by assuming that the appropriate singular field is a 4D Hadamard Green’s function 
on the constant time Riemannian slice. It is well known that the Hadamard representation of 
a Green’s function involves an arbitrary regular biscalar Wq{x,x'), whose coincidence limit w(x) 
appears in the expression for the self-force. We develop an axiomatic approach to reduce this 
arbitrary function to a single arbitrary dimensionless coefficient. We show that in the context of 
this approach to regularization, the self-force does not depend on any undetermined length-scale 
and need not depend on the internal structure of the charge. 


I. INTRODUCTION 

The problem of motion in General Relativity is ex¬ 
tremely complicated owing to the non-linearity of Ein¬ 
stein’s equations. However, when there is a natural sepa¬ 
ration of scales, for example an extreme mass ratio binary 
black hole system, the problem is amenable to a pertur¬ 
bative analysis. To zeroth order in such a perturbative 
scheme, the small body is described by a point particle 
moving on a geodesic. Going beyond the geodesic ap¬ 
proximation by including higher order corrections in the 
perturbative analysis can be interpreted as arising from 
the self-interaction of the point particle with its own field; 
this is the so-called self-force problem. In four dimen¬ 
sions, the self-force has enjoyed a long and fruitful history 
and is very well understood both at a formal level and 
at a practical computational level (see Ref. [T] for a com¬ 
prehensive review of the self-force problem and Ref. [2] 
for a review of computational approaches). Much of the 
attention in the literature is devoted to the development 
of a regularization prescription for curing the divergences 
that inevitably arise when working in the point-particle 
limit. Curing this kind of pathology is standard fare in 
quantum field theory where the infinities are absorbed 
into a renormalization of the constants which are then 
determined by observation. The diSiculty in the classical 
theory is constructing the unique, finite self-force that 
corresponds to the self-force on a finite-size body in the 
point-particle limit. There are a number of equivalent 
ways to regularize the self-force in curved spacetime dat¬ 
ing back to the seminal paper of DeWitt and Brehme [3] 
who computed the formal expression for the regularized 
electromagnetic self-force, which was later corrected by 
Hobbs [4]. Their approach relied heavily on the covari¬ 
ant decomposition of the Green’s function for the elec¬ 
tromagnetic vector potential into a “direct” and “tail” 
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piece. The direct part, which they take to be the aver¬ 
age of the advanced and retarded field, is singular and 
has support only on the light-cone while the tail part is 
regular with support only inside the light-cone. 

The most transparent derivation of the self-force is 
furnished by the method of matched asymptotic expan¬ 
sions which was utilized by Mino, Sasaki and Tanaka 
to formally derive the self-force for a point mass in 
curved spacetime (see also Refs. mm- A complemen¬ 
tary derivation of the gravitational self-force was offered 
by Quinn and Wald [7] which was based on a simple set 
of physically motivated axioms. Hence, the equations 
of motion including first-order self-force effects are com¬ 
monly referred to as the “MiSaTaQuWa” equations. In 
the same paper, Quinn and Wald also derive the electro¬ 
magnetic self-force and their result agreed with that of 
DeWitt and Brehme [3] (taking into account the correc¬ 
tion by Hobbs i)- This simple axiomatic approach was 
later extended to the scalar self-force by Quinn [S] . Each 
of the aforementioned calculations adopted the regular¬ 
ization scheme of DeWitt and Brehme [3] . An alternative 
regularization prescription was offered by Detweiler and 
Whiting [S] who instead decomposed the physical field 
into a “singular” and “regular” field and showed that 
the regular field is a solution to the homogeneous (source- 
free) wave equation and hence is smooth at the particle’s 
location. This decomposition leads to a more natural 
interpretation of the self-force in terms of a particle’s in¬ 
teraction with an external source-free field. Moreover, 
an axiomatic construction of the Detweiler-Whiting sin¬ 
gular field was developed in Ref. [T] by demanding that 
the singular Green’s function is a symmetric, inhomoge¬ 
neous solution to the wave equation with zero support 
inside the light-cone. 

The gravitational self-force was given its most rigorous 
treatment by Gralla and Wald [10] who considered a one- 
parameter family of metrics describing an extended body 
that is scaled down to zero size and mass in an asymp¬ 
totically self-similar manner. This derivation is essen¬ 
tially a more rigorous version of the matched asymptotic 
derivation given by Mino, Sasaki and Tanaka [^ . A rigor- 
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ous derivation of the electromagnetic self-force was later 
given by Gralla, Harte and Wald m- 

The self-force itself is not directly observable; instead 
the quantity of interest is the gravitational waveform pro¬ 
duced by the emitted gravitational waves, for example 
by a particle inspiralling into a black hole. Computing 
the waveform itself requires solving the coupled problem 
of particle motion and radiation generation. If one uses 
straightforward black hole perturbation theory, geodesic 
particle motion is associated with (i.e. arises at the same 
order as) first order metric perturbations. Accelerated 
particle motion and the hrst order self force are associ¬ 
ated with second order metric perturbations. This per¬ 
turbation theory can self-consistently describe acceler¬ 
ated particle motion and the associated waveforms that 
are generated, but only over times short compared to the 
dephasing time (the timescale over which a geodesic or¬ 
bit and the true orbit become out of phase by ^ 1 cycle). 
Thus, computing waveforms that are accurate over a ra¬ 
diation reaction timescale requires going beyond this kind 
of perturbation theory. Two complementary approaches 
for achieving this are the self-consistent worldline ap¬ 
proach of Pound [12] and the two-timescale expansion 
technique of Hinderer and Flanagan m- 

Despite the comprehensive work and progress in the 
self-force problem in four dimensions, the situation re¬ 
mains underdeveloped in higher dimensions. Obviously 
the intense endeavor over the last few decades has been 
driven by modeling extreme-mass-ratio inspirals and 
hence the self-force in a more abstract context has re¬ 
ceived little attention. Moreover, there is as yet no 
known general prescription or expression for self-forces 
in higher dimensions. Despite the lack of formal under¬ 
pinning, there has been considerable recent interest in 
the self-force in higher dimensions [HHiij yielding some 
very unexpected results in odd dimensions. In particu¬ 
lar, Beach, Poisson and Nickel [14] calculated the self¬ 
force on a static scalar and electric charge in a 5D black 
hole spacetime and found that the result depended on 
the radius of a sphere centered at the charge, which they 
interpret as the radius of the particle. Moreover, Frolov 
and Zelnikov [TS] calculated the self-force on a uniformly 
accelerating charge in Minkowski spacetime (or equiva¬ 
lently a static charge in Rindler spacetime) and their re¬ 
sult depends on an undetermined infra-red cutoff length- 
scale which the authors postulate could be the scale over 
which the homogeneous gravitational field approximation 
is valid. 

In this paper, we will revisit the calculation of a static 
electric and scalar charge in a 5D black hole spacetime. 
We derive closed-form representations of the electrostatic 
field and the static scalar field. In the absence of a rigor¬ 
ous derivation of the self-force in five dimensions, we at¬ 
tempt to use the Hadamard regularization prescription to 
compute a locally constructed singular field to subtract 
from the self-field. For a static charge in a static 5D 
spacetime, the held equations reduce to an elliptic wave 
equation on a 4D Riemannian manifold. It is well-known 


that the Hadamard form of the Green’s function for this 
wave equation involves an arbitrary bi-scalar Wq{x,x') 
that is undetermined by the local Hadamard expansion. 
For the physical held, this arbitrariness is a failure of 
a local expansion to encode global information such as 
boundary conditions or information about the quantum 
state in the quantum theory. For the singular held, we 
develop a set of axioms for constraining the coincidence 
limit of this bi-scalar, valid for static charges in an ar¬ 
bitrary static 5D spacetime. Remarkably, our axioms 
reduce a functions worth of arbitrariness down to a sin¬ 
gle arbitrary dimensionless coefficient. Using our closed- 
form expression for the physical held and our axiomatic 
construction of the singular held allows us to calculate 
simple closed-form expressions for the self-force in terms 
of this arbitrary coefficient. Our regularization scheme 
yields a self-force that does not depend on any undeter¬ 
mined length-scale such as an infra-red cutoff, unlike the 
calculation of Ref. [T^ , nor does it depend on some par¬ 
ticular model for the point particle, contrary to the result 
of Ref. [II]. 

II. THE ELECTROSTATIC FIELD 
A. The Wave Equation 

The (n-|-2)-dimensional analogue of the Schwarzschild 
black hole, known as the Schwarzschild-Tangherlini met¬ 
ric [TH], is given by 

ds^ = -f{r)dt^ + j^dr'^ + r^dfl^, (2.1) 

where /(r) = 1 — (rH/r)"”^ with the horizon radius 
and is the line element on a unit n-sphere which may 
be dehned inductively by 

dfl^ = VlABdx^dx^ = d9i + sin^ (2.2) 

The angular variables are x^ = {0i,..., 9n-i, (/>} with 9 G 
[0,7r) and (/> G [0,27r). We can re-write the metric in 
terms of a lapse function and a Euclidean metric on a 
constant time hypersurface by 

ds^ = gabdx°‘dx^ = —N'^dt^ -\- hapdx°‘dx^, (2.3) 

where Latin indices a,b — 0, ..,n + 1 are used for space- 
time components and Greek indices a, fd = l,...,n-|-l 
refer to tensor components on the constant time Rie¬ 
mannian slice. We are interested in computing the self¬ 
force on a static electric charge in this spacetime. The 
only non-trivial component of the vector potential is 
$ := and Maxwell’s equations reduce to the (n-|-l)- 
dimensional Helmholtz equation on the Euclidean metric 

^a/3 

{v2-A“a<,}$(a^) = D„/j‘ (2.4) 
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where 

1 27r(”+^^/^ 

~ p^ rt+l ^ ~ 

and V^ is the Laplace operator on the metric The 
charge density for the point charge e at position Xq is 
given by the Dirac delta distribution 

j* =-^J(r - ro)<5(D,f2o), (2-6) 


Since the spacetime is spherically symmetric, the ra¬ 
dial Green’s function gi{r, r') is independent of the mode 
number j and hence we can apply the addition theorem 

A(Z,n)-l 

^z(cos7n) = -^ E (2.13) 

where 7 „ is the geodesic distance between two points on 
S" and Vi is the generalized Legendre polynomial which 
is a solution to the differential equation 


where (5(D, Dg) = S^(x'^ — x^) / ^J\^Iab\ is the invari¬ 
ant Dirac delta distribution on an n-sphere. For a static 
particle in a static, spherically, symmetric spacetime the 
only non-vanishing component of the self-force is in the 
radial direction and is formally given by 

F'' = e^odMxo), (2.7) 

where /o = /(r’o). This expression is meaningless as 
it stands since the field diverges at the position of the 
charge owing to the distributional nature of the point- 
particle source. Hence the gradient of the field must 
be regularized before evaluating at the position of the 
charge, which we discuss further in Sec. |IV A| 

Now, for the particular spacetime under consideration, 
the Green’s function for the wave equation is 

{ V2 - A‘^d^}G(x, x') = _ r')5{n, D'), 

( 2 . 8 ) 


which can be recast as 


+ G{x.x') = 


-D, 


5{r - r')5{n,VL') 
^/n-2yjpy 
(2.9) 


■ (P d 

+ cot 7 „-- \-l{l + n 

La7„ d'jn 

normalized so that 


1) 'P;(C0S7„) 


0, (2.14) 


Vi (cos 7 „)T’/' (cos 7 „) sin” 7 „ djn 


Hn—1 


6iu 

A(Z, n) 
(2.15) 


The electrostatic Green’s function now assumes the form 

OO 

G{x,x') ='^A{l,n)Vi{cos-fn)gi{r,r'), (2.16) 
z=o 


which upon substitution into the wave equation gives the 
following inhomogeneous ordinary differential equation 
satisfied by gi{r, r'), 


r —— -I- nr-- 

ar^ dr 


d l{l + n — 1) 


/ 


gi{r, r') 


5{r — r') 


(2.17) 


Rewriting this equation in terms of the dimensionless ra¬ 
dius 


C = 2(r/rH)”-i - 1, (2.18) 


with the Laplace operator on the n-sphere whose 
eigenfunctions are the generalized spherical harmonics 
m with eigenvalues l{l + n — 1). For a given 

I, there are 


A{1, n) 


(21 + n — 1) {I + n — 2)! 
(n- 1)! h 


( 2 . 10 ) 


linearly independent harmonics and the angular Dirac 
delta distribution may be decomposed in this basis as 


leads to 



Z(/ -I- n 
{n- 


2 <5(e-o 


(«-!)<■' y7(F) ■ 


(2.19) 


Following standard techniques (see [5D] for example), we 
can construct a solution to this inhomogeneous equation 
as a normalized product of homogeneous solutions. 


OO A(i,n) —1 

d(H,D') = E E y/j(^)yij(^/)^ (2.11) 

1=0 j=0 

This suggests the following mode decomposition for the 
Green’s function 

cxo A(l,n) — 1 

G(x,x') = n^J^ V‘’Pn)Y‘’PQ')gi(r,r'). 

1=0 j=0 

( 2 . 12 ) 




2 ^r^(g<)^pHe>) 

(n-l)rr'y7(F) Gi 

( 2 . 20 ) 


where and 'I>[^^(^) are linearly independent so¬ 

lutions to the homogeneous equation satisfying regu¬ 
larity on the horizon and at infinity, respectively, and 
Gi is determined by the Wronskian of these solutions. 
We have adopted the notation = min{^,^'} and 
= max{^,^'}. The linearly independent solutions 
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of the homogeneous equation are •\/ 7 ^//(^„_i )(0 
V 7 Qi^/(„_i)(C) where Pjfiz) and Q'iiz) are the associ¬ 
ated Legendre functions of the first and second kind, 
respectively. The former is the solution regular on the 
horizon while the latter is regular at infinity, i.e., 

( 2 . 21 ) 

where A = l/{n — 1). 

However, for I = 0, these regularity conditions do not 
fix the solutions since both are everywhere regular and 
the choice is determined by enforcing the total charge as 


measured by an observer at infinity to be the charge of 
the point particle. This ambiguity is well known pTl [22] 
and the appropriate choice of solutions is a constant for 
the inner solution and l/(^-|-l) for the outer solution. We 
prefer to express these in terms of associated Legendre 
functions, i.e., 

= ^/fQoiO- ( 2 . 22 ) 

Noting that the normalization constant is simply Ci = 
— 1 for all I, our mode-sum representation of the Green’s 
function is 


G{x, x') = —~ 


2^7 


(n-l)rS 


■E 

1=0 


A{l,n)Vi{cosjn) Px \U)Qii^>) + 


4V/ 


(n - l)rl 


xtQUOQUO- 


(2.23) 


r 


Finally, the convolution of the Green’s function with the 
source term appearing in Eq. (|2.4|) yields 


(2.24) 


4>(a;) = -ex/foG(x,Xo), 


which results in the mode-sum representation of the elec¬ 
trostatic potential 

®(®) = '^Hl,n)ri{cosjn)Px^{U)Qi{i>) 


{n — l)r 


1=0 


4e 


(n-l)rr' (^ + 1)(6 + 1 ) 


(2.25) 


where we have used the fact that / = (C — 1 )/(C + 1 ) 
and Qo(C) = —— 1. In the four-dimensional 
Schwarzschild spacetime, we can sum these modes to re¬ 
trieve the Gopson-Linet [2U [53] closed-form electrostatic 
field. 


B. Closed Form Static Field in Five Dimensions 

In this section, we specialize to the case of five dimen¬ 
sions (n = 3). We will derive the closed-form representa¬ 
tion of the electrostatic field by summing the mode-sum 
representation derived in the previous section. 

For n = 3, the generalized Legendre function has the 


simple trigonometric representation 

sin(; -I- 1)73 


T"/(cos 73 ) = 


{I + 1 ) sin 73 ’ 


(2.26) 


which is easily verified by substituting into Eq. (2.14|. 


Henceforth, we shall drop the subscript on 73 for typo¬ 
graphical convenience. The mode-sum representation of 
the field is now given by 

+ l)7.P/72(C<)Qi%(?>) 


sm 7 


2r'^rQ 


1=0 


(2.27) 


Now consider the standard Legendre addition formula 


Q.(eeo-(C"-I)'/'(?o'-l)'/'cos4-) 

00 

= ^ e-^'^-p-\U)Ql{i>)cosk^ (2.28) 

k— — oo 

valid for all 1 / ^ —1, —2, —3,..., and taking the Fourier 
inverse allows us to express a product of associated Leg¬ 
endre functions in terms of a single Legendre function. 


e-^'^^p-H^<)Qti^>) = cosfc4/Q,(eeo - (.e - I)'^'(eo " I)'/'cos 4/)d4/. 


(2.29) 
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Employing this identity in our expression for the field above yields 

eVfvTo 


= — 


27rr| sin 7 Jq 


/ cos^'^(/ + 1 ) sin(Z + 1 ) 7 Q// 2 (CCo - “ 1 )^^^ cos - ^ 7 ^. 

•^0 1=0 '’0 


In Ref. [25], we find the summation formula 


E ,, , „ , , i /TT 2 cos 7 \ 

■=»(/ + 1)7 Q,M-> = (2,_2cos27)-/d 2 + (2z - 2cos27)-/d ' 


Differentiating with respect to 7 gives 


+ 1 ) sin(; + 1 ) 7 Qz/ 2 ( 2 ) = 


sin 7 


+ 


TT Sin 7 cos 7 


^ sin7COS7arctan( ^J_ 77^7 


;=o 


(z-cos 27) v^(z - COS 27)3/2 


(z — COS 27 ) 3/2 


r 


(2.30) 


(2.31) 


(2.32) 


which permits us to sum the modes in our electrostatic 
potential to obtain the integral representation 


<i)(a;) = — 


’■VfVfi 


2tt r 
■\/ 2 cos 7 

o3/2 


(A r ( 

I Jo V 


1 ^ COS 7 


arctan 


P P' 

1 —b 

V 2 cos 7 / / 


,3/2 


cos 'I' d'h — 


2 r 2 r 


(2.33) 


where 


p = - cos 27 - (e^ - 1)1/2 (^2 _ 1)1/2^03^^ ( 2 . 34 ) 


In arriving at (2.33), we have restricted attention to 0 < 

(2.35) 


7 < 7r/2 and made use of the identity, 

arctan(a;) = f — arctan(l/a;), a; > 0 . 
Defining 

_ iio - cos 27 
X- (^2_1)1/2(^2 _i)1 /2> 

we can recast our integrals as 

e v7V7o ( 1 


(2.36) 


<I>(a:) = — 


27rr2 I (^2 _ l)l/2(^2 _ 1)1/2 


X 


e) 




TT cos 7 


(^2 - l)3/4(Co2 - 1)3/4 


X 


e) 


72 


COS 7 


(7 - l)3/4(e§ - 1)3/4 


X 


(e) 


_e7 

2r^r 


(2.37) 


where 


x7 = 




cos 


b (x-cos^-) 

p27r 




COS ^ 


fo 


(x — cos 




7-(®) _ 

-^3 ~ 


L 


2-77 arctan 


V 2 cos 7 ^ 


cos'I' 


(X — cos 5')3/2 


The hrst two of these are straight-forward to integrate; 
introducing the notation 


3± 


= x± 7 x^ - 1 , 


(2.39) 


then we have 


i 7 = 2 . 


x 7 ^ = -2 


77^’ 

d /■ 2 ’^ 


COS T 




dxJo (x-cos«')i/2 

dx V 7 ^ 

where E{z), K{z) are the complete Elliptic integral func¬ 
tions of the first and second kind, respectively. In arriv¬ 
ing at the last line above, we made use of standard iden¬ 
tities for derivatives of Elliptic integrals [24] and we also 
used the fact that z_ = l/z+. 

The third of these integrals seems troublesome at first 
glance, but can be done by moving to the complex plane 
and treating the branch points with sufficient care. We 
make the following transformation 




dz 

d"^ =— and cos'^ = ^{z1/z). 
iz ^ 


(2.41) 

In terms of z, the integral goes around the unit circle in 
the complex plane and may be written as 


= *72 


■d^'. 


/(7(7 + i) 

(z-Z_) 2 (z-Z+> 


rdz, 


(2.42) 


(2.38) where 
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/W 





2 cos 7 V z 


(2.43) 


Im(z) 



FIG. 1: Plot of the branch cuts and the contour of integration 
iorT^^\ 


and z± are defined in Eq. (2.39). One can show that f{z) 
is holomorphic in a neighborhood of z± but has branch 
points at z = 0 and at 


z± = X ± Vx^ - 1 


where 


^^0 + 1 


(2.44) 


(2.45) 


We note that z_ < and z_|_ > and hence introduc¬ 
ing branch cuts on the real axis from [ 0 , z_] and [z+, oo) 
yields a single-valued holomorphic integrand except at 
isolated simple poles at z±. Of these poles, only 2 ;_ lies 
inside the unit circle. 

Now we consider the integral around the deformed con¬ 
tour r which consists of (see Fig. [^: 

(i) a circular arc Fi traced counterclockwise from 

— Vl — — *e to —-s/l — -I- ie, 

(ii) a horizontal line segment Ii from — \/l — 4- ie to 

z- -I- ie, 


(iii) a semi-circle Fj of radius e centered at and traced 
clockwise and 


(iv) a horizontal line segment I 2 from z_ — ie to 

— Vl — — *£■ 

Hence the integral around the contour F is given schemat¬ 
ically by 


= + 


-f 


’h 


'h 


which can be evaluated by the Cauchy Residue Theorem 
yielding 


= -TTl- 


_.(^2_i)1/4(^2_ 1)1/4 


X 


cos 7 


Vx^ -1 


(2.47) 


In the limit as e —)■ 0 the integral around F 1 tends to the 
closed contour around the unit circle that we required 
and so we obtain 


= lim 

£->■0 


= —7^^- 


.(^2_i)1/4(^2_i)1/4 ^ 


\/x^ -1 


cos 7 

— lim 
£->■0 

(2.48) 

Following standard text-book techniques, it is straight¬ 
forward to show that the integral around F^ vanishes as 
e —?► 0. Moreover, in this limit the line integrals Ii and I 2 
cancel in the left-half complex plane since the integrand 
is single-valued there, but we pick up a contribution from 
these line integrals across the cut since the arctanh is dis¬ 
continuous across the cut and the difference above and 
below is ITT. Therefore, we have 


lim ( / + [ ) = lim ( [ 


Z— -\-ie 


+ 


0+ie 

2 _ 


= —ZTT 


^ 0 —ie 


/ Z- —le 

+ 1 


yjx{x^ — 2xx + 1 )^/^ 


dx. 

(2.49) 

This integral can be written in terms of elliptic integrals 
as 

f x^ + 1 _ _ (xEiib z ) 

Jo V^ix^-2xx + ir/^ (x^-l)V^ 

-Vx^ - lE(V’,2r-)) 

I (^^ - - 1 )^/^ x(x - Z-) 

X 2-1 ’ 

(2.50) 


cos 7 


where 


Ip = arcsin(-\/z_ z+). 


(2.51) 


Combining Eqs. (|2.48|)-(|2.50|) with Eq. (|2.42|) yields 
7r-\/2 


t(®) _ 

Xo — 


x'-l 


2 v^(x£’(V'. Z-) - 'Jx^ - 1 ^(V', 2 ;-)) 

(^2_i)1/4(^2_i)1/4 


cos 7 


x(z_ - Z-) 

(2.52) 


(2.46) 









































7 


Finally, the electrostatic field is given by 

2 x^ - X 5- - 1 


<i)(a;) = — 


+ 1)(6 + i)^x^ -1 
2,/ircc»7 ,2{Kiz.)- 


'/x^ -1 


X 


(e - - 1 )'/" 


\/x^ -1 


- F{^, z.) - 


X 


\Jx^ - 1 




£;(z_) 


(2.53) 


where 7, CiXi Zj-, Zj- and ip are defined in Eqs. (2.13), 


(2.18), (2.36), (2.39), (2.44) and (2.51), respectively. 


Now, recall that, since we have employed the identity 


(|2.35|) to compute this expression is valid only for 


0 < 7 < 7 r/ 2 . We can analytically continue the Elliptic 
integral functions that appear in our expression for 
which results in the following representation for the held, 
valid over the entire range 0 < 7 < tt. 


<i)(a;) = — 


r2(^ + l)(^o + l)\/x^ - 1 

2 COS 7 


2 x^ - X g_ - 1 

Vx^ - 1 


—====£(7_)j - (2e(cos7) - 7_) 


X 


y/x^ - 1 


E{i,, z.) 


(2.54) 


where Q{z) is the Heaviside step function. We have 
checked that this held does indeed satisfy the wave equa¬ 
tion and that it agrees numerically with its corresponding 
mode-sum representation. 


III. THE SCALAR FIELD 

In this section, we apply the previous analysis to a 
static scalar charge. The calculation is much the same as 
the electrostatic case and so we proceed with less details 
than before. 


A. The Wave Equation 

A static scalar particle with charge q at position Xq in 
a static, spherically symmetric (n-|- 2 )-dimensional space- 
time satishes the (n-l-l)-dimensional Helmholtz equation 


where 


{V^ -I- A°‘da}‘p{x) = -Xln n 

^i = q - ro)S{n, Hq). 

'^n 


(3.1) 


The vector held appearing in the potential is the same 
as in the electrostatic case (2.5) but note that the poten¬ 


tial has the opposite sign. The scalar self-force arising 
from the static scalar charge interacting with its own 
held has only a component in the radial direction and 
is formally given by 


F'" = qfodr(p{xo). 


(3.3) 


As before, we need to regularize the gradient of the held 
before evaluating at the location of the charge. The cor¬ 
responding Green’s function for the static scalar wave 
equation is 


= -a 


{V2+A“5„}G(a:,a:') 

,^/®^(r-r')<5(H,H') 


(3.4) 


which is precisely Eq. ( |2.8[ ) with the transformation 
A“ —>■ — A“. For the Schwarzschild-Tangherlini space- 
time, this yields 


( 2 ( 

I- 


f ) dr ^ f " 


|G(a;,a;') 


Y^/rL — 2 


\/J¥) 


S{r — r')S{n,Xl'). (3.5) 


The Green’s function is again decomposed in a basis of 
generalized spherical harmonics as in Eq. (2.16), where 
now the radial part gi{r,r') satishes 


f 2 V ( . rfo\ d l{l + n-l)^ , ,, 

F ^+d"+Tl*- -f — 


5{r — r') 


(3.6) 


In terms of the dimensionless radius ^ dehned in 
Eq. (|2.18|), this is 


l{l + n — 1 ) 
. V''’ ^ d^J (n — 1)2 






2 vW 0 


(n — 1)7 


^<5(^-r). (3.7) 


The solutions of the corresponding homogeneous equa¬ 
tion are Legendre functions of degree A = Ijin — 1). We 
choose the Legendre function of the hrst kind P\{^) to 
be the inner solution regular at the horizon and the Leg¬ 
endre function of the second kind Q\{^) to be the outer 
solution regular at inhnity. Unlike the electrostatic case, 
there is no monopole ambiguity and so the mode-sum 
representation of the Green’s function is 

G{x,x') = 7-^-™i^^A(Fn)Xz(cos7„)PA(C<)QA(C>)- 


(n - l)rS- 


1=0 


(3.2) 


(3.8) 

















































To obtain the field at x due to a static charge at a;o, 
we convolve the Green’s function with the scalar source 
which results in 

9 /~f~ ^ 

v{x) = - - ^ °_1 ^A(^,n)iP;(cOS7„)PA(^<)QA(C>)- 

(3.9) 


B. Closed Form Static Field in Five Dimensions 


sum representation of the static scalar field 


>f{x) = 2 ^ + sin(Z + l) 7 ^)/ 2 (C<)< 9 // 2 (C>)- 


sin 7 




(3.10) 


For n = 3, we adopt the representation (2.26) for the 


generalized Legendre function which leads to the mode- From Eq. (2.29), we have that 


1 

Pl/2{U)Ql/2{^>) = QlM^^O - ie - l)'/"(eo - l)'/"cOSVl/)dvI/, 


which allows us to recast the mode-sum as 

/‘27r oo 

9 V/o 


rf^ pZTT 7 ^ 

‘p(x) = ^ ^ 2 ° / ++- (^^ - l)^/^(^o - l)^^^cos^')dT. 

Sin 7 Jo ^^0 


(3.11) 


(3.12) 


r 


As in the electrostatic case, the sum can be performed 
using Eq. (2.32) which results in the integral representa¬ 
tion 

/ N (iVTo f ^ C 0 S 7 


StitI 7o p3/2 

cos 7 arctan( 7 ^) 


.3/2 




(3.13) 


where p is defined by Eq. (2.34). We write these integrals 
in terms of x El'S 


(p{x) = 


qVfo ( 1 

27 rri 172-1)1/2(^2-01/2-1 


■T 


(s) 


7 r -\/2 cos 7 


(e - - 1)^/4 


t'-' 

Xo 


V2 


COS 7 


ie - i)3/4(^o" -1)3/" 


Xq 


(3.14) 


where 

x(^) = r _^_ 

^ Jo ix-COS'S) 




(x — cos 4>)3/2 


d^'. 


i7 = 


27r arctan[l ^— "7 —(x ~ cosT)^/^] 


x/2 cos 7 

(x — cos 'll)3/2 


-d^'. 

(3.15) 


The first integral can be performed in terms of elemen¬ 
tary functions while the second results in a combination 
of complete Elliptic integral functions. Specifically, we 
have 


tC®)_ 

-^1 “ 


27r 


7x^-1 


= -4v/2 


Vx^ - 1 


(iX(z_) 




7x^-1 


E(z_) 


(3.16) 


where z± are given by Eq. (2.39). The remaining integral 


(s') . 

I 3 is done by taking the integral around the unit circle 
in the complex plane 


= * 2^2 / 


■■fiz) 


I {z-z+)^{z-z-} 


•dz, 


(3.17) 


where f{z) is defined in Eq. (2.43). We consider the same 
deformed contour as in the electrostatic case and again 
we pick up a contribution only from the pole at z = z_ 
and a contribution across the branch cut which lies along 
[0, z_]. The result is 


q-{s) _ TTV^ 

3 “X "-1 


2v^(z+E(V’, Z-) - 7x^ - z-)'^ 

(7-1)1/4(^2-1)1/4 


COS 7 


-(i_-z_) 


(3.18) 


where ip is given by Eq. (2.51). 






































9 


Finally, the static scalar field is given by 


(p{x) = 


<lVfo 




2 cosy^^ 


(e - 


(^F{y, Z-) 


Z+ 


\/x^ - 1 


!(iF(z_) 


E{ip, z_) 


X- z- 


\/x^ - 1 




(3.19) 


where 7 , % , Z j-, Zj- a nd -0 are de fined in Eqs. (2.13), 

(2.18), ( |2.36 ), (2.39), (2.44) and (2.51), respectively. 
Again, this expression is valid only for 0 < 7 < 7 r /2 
but can be analytically continued as before, resulting in 


higher dimensional analogue of the Detweiler-Whiting 
Green’s function, namely, 

Gow(a:,a;') = | (^)d/2 

(4.2) 

This Green’s function can be constructed axiomatically 
[T] by demanding that the appropriate parametrix satisfy 

1 . (□ - C R)Gsix, x') = -nd- 2 S{x, x'), 

2. Gsix,x') = Gs{x',x), 


ip{x) = 




rue-iy/yeo-^y/^Vx^ 


X- z- 

\lx^ - 1 


2 cos 7 ^fzZ \ , s f , . 

(P-i)i/4(|_i)i/4 j^e(°”r)(A'(^-) 

r)^ E(=-)) - (2e(cos7) - l)(f 7 ’, 2-) 




\/x^ -1 


E{y, z-) 


(3.20) 


We have checked that this field satisfies the static 
scalar wave equation and agrees numerically with its cor¬ 
responding mode-sum representation. 


IV. THE SINGULAR FIELD FOR STATIC 
CHARGES 


A. Hadamard Green’s Functions 

The formal expression for the self-force on a point 
charge involves taking the gradient of a field evaluated 
at the location of the charge. This is clearly divergent 
as implied by the delta distribution source and therefore 
requires regularization. In four dimensions, the most ele¬ 
gant regularization prescription that results in the correct 
self-force is the Detweiler-Whiting [S] construction of the 
singular field which yields, upon subtraction from the re¬ 
tarded field, a regular homogeneous solution to the wave 
equation. The Detweiler-Whiting Green’s function is 

G^^{x,x') = ^(^A^/2(a;,a;')^(cr) - 0(cr) F(x, a;')) . 

(4.1) 

where a{x,x') is Synge’s |26j world function which is 
half the square of the geodesic distance between x and 
x', 0(z) is the step function, A(a:,a;') is the Van Vleck- 
Morette determinant and V(x, x') is a regular, symmetric 
bi-solution of the homogeneous equation. In even dimen¬ 
sions, it is a straight-forward matter to construct the 


3. Gs(a;, x') = 0, if a{x,x') <Q. 

In odd dimensions, by contrast, it seems likely that no 
Green’s function satisfying these properties exists, even 
in Minkowski spacetime; a variety of simple choices of 
ansatz can be shown to not satisfy the properties above. 

If we are only interested in a static charge in a static 
spacetime, then the wave equation reduces to an ellip¬ 
tic equation on a Riemannian manifold, which has a 
unique solution subject to boundary conditions. Al¬ 
though the Green’s function for the physical field is sen¬ 
sitive to global properties such as boundary conditions, 
within some local neighborhood all Green’s functions for 
all choices of boundary conditions are described by the 
same universal Hadamard form, which we define below. 
Each propagator within this family of Hadamard Green’s 
functions is parametrized by a particular choice for some 
regular biscalar W^x^x'). It is the difference between 
this biscalar Wq{x, x') for the physical field and the singu¬ 
lar field that gives rise to the self-force. We assume that 
our singular field constructed from a Hadamard Green’s 
function (for some judicious choice of Wo{x, x')) does not 
exert a force on the charge and hence subtracting from 
the physical field and taking the gradient leads to the cor¬ 
rect self-force. This assumption can be partially justified 
on the basis of Harte’s formalism [271 - 130] for self-fields on 
extended bodies wherein it is proved that the self-force 
contribution from a singular field obtained from a geo¬ 
metrically constructed symmetric two-point function acts 
only to renormalize the moments of the body. Hence our 
results for the self-force ought to be correct up to such 
renormalizations. 

In his seminal lectures, Hadamard m constructed a 
local solution to an arbitrary second order linear partial 
differential equation by assuming the solution has a se¬ 
ries expansion in Synge’s world function (j{x,x'). The 
expansion is only valid in a neighborhood where x and x' 
are connected by a unique geodesic, the so-called normal 
convex neighborhood. We present only a brief overview 
here and we refer the reader to the detailed description of 
Hadamard renormalization in arbitrary dimensions given 
in Ref. |5^ . 
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In d odd dimensions, the universal Hadamard form for 
a Green’s function takes the form 

Gh{x,x) ^ +W{x,x), (4.3) 

where U{x, x') and W (x, x') are regular biscalars and U 
satishes coincidence boundary condition [U] = 1 (we have 
adopted square brackets around a bi-tensor to denote the 
coincidence limit x —)■ x'). Substituting the ansatz 

OO 

U{x,x') = ^ Uk{x,x')a^ 

k^O 

OO 

W{x, x') = ^k{x, x')a'^ (4.4) 

k=0 

into the wave equation yields a set of recursive trans¬ 
port equations for the coefhcients and 144- The re¬ 
cursion relations completely determine all of the Uk co¬ 
efficients and also determine all of the Wk coefficients 
except for Wq which remains arbitrary. This implies 
that the Hadamard representation of a Green’s function 
is not completely determined by a local expansion. For 
the physical field, Wq encodes global information such as 
boundary conditions. In quantum held theory, for the 
Hadamard representation of the two-point function, Wq 
encodes the information about the quantum state of the 
system. 

In d even dimensions, the Hadamard form is given by 

Gh( x,x') -- ^^d/ 2-1 +y{x,x')\oga + W{x,x') (4.5) 

where U{x,x'), V(x,x') and IF(x,x') are regular bis¬ 
calars. One can make the log term explicitly dimen¬ 
sionless by introducing a length scale s, and writing it 
as V\oga/s^ but this length scale can be absorbed into 
a redefinition of IF(x,x')- The biscalar U{x,x') again 
satisfies the boundary condition \U] = I, but unlike the 
odd-dimensional case now possesses a finite series expan¬ 
sion 

d/ 2-2 

U(x,x')= Uk[x,x')u^. (4-6) 

k=0 

Assuming the series expansions 

OO 

H(x,x') = ^ Vk{x,x')(j^, 

OO 

W{x,x')=Y,Wu{x,x')a\ (4.7) 

fc =0 

and substituting into the wave equation determines the 
Uk and 14 completely. As before, this determines all of 
the Wk coefficients except for Wq. 

If the wave equation is self-adjoint, say, 

- P(x)|Gh(x,x') = -Od -2 S{x,x') (4.8) 


where P{x) is a potential that does not contain any 
derivative operators, then t7(x,x'), V{x,x') and lT(x,x') 
are symmetric biscalars. The biscalar IF(x, x') is then 
a bi-solution of the homogeneous wave equation in odd 
dimensions while V{x,x') is a bi-solution of the homoge¬ 
neous wave equation in even dimensions. Furthermore, 
we have that (see, for example, [331 [Mj) 

Uq{x,x') = = 1-1- ^Ra'i3'cr’°' 

cr’^ 0 -’^ -f o(cr^) (4.9) 

while Vq and Wq possess covariant Taylor expansions of 
the form 

Vo(x,x') = v(x') - ^V;a'(x')a’°' + 0(a), 

Wo(x,x') = w(x') - ^W;a'(x')a'‘^ -I-O(ct), (4.10) 

where w(x') = [y(x,x')] and w(x') = [IF(x,x')]- Explicit 
expressions for v(x') are obtained by taking coincidence 
limits of its transport equation, for example, for d = 4 
we have v(x') = ^(P(x') — ^R(x')). 

B. Constraining w(x) for a Static Charge in a 5D 
Static Spacetime 

For a static charge located at Xg in a 4D static space- 
time, it has been shown |3S] that the Detweiler-Whiting 
Green’s function is equivalent to the three-dimensional 
Hadamard Green’s function on a constant time slice with 
w{xo) = 0, at least up to the order required for the com¬ 
putation of the self-force (assuming the singular field is 
derived from a symmetric Hadamard Green’s function, 
then only w(xq) and its derivative appear in the self¬ 
force). For ultrastatic spacetimes (gtt = —1), they are 
equivalent up to all orders with W(x,Xo) = 0. In these 
cases, it is possible to determine w(x) by direct com¬ 
parison with the known Detweiler-Whiting field. In this 
section, we address whether we can determine w(x) with¬ 
out knowing the corresponding Detweiler-Whiting field. 
We note that in Ref. the authors adopt a defini¬ 
tion for the self-force that involves a spherical averaging 
procedure and the smooth part of the Green’s function 
does not contribute to the self-force after this averaging 
is performed. Thus, their analysis effectively corresponds 
to choosing the singular field for which Ho(a;,a:o) = 0. 
However, such a spherical averaging results in a self-force 
that depends on the details of the model charge. For ex¬ 
ample, one can model the point particle as a sphere of 
constant proper radius centered at the charge’s location 
or alternatively as a sphere of constant proper volume 
centered at the charge’s location; the resultant self-force 
is different for each model. Moreover, the self-force di¬ 
verges as the radius of this sphere shrinks to zero. 

Gonsider now a static electric charge in a d-dimensional 
static spacetime. We decompose the metric into a lapse 
function N and a (d — l)-dimensional Riemannian met¬ 
ric hap. The only non-trivial component of the vector 
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potential is the time component <i)t = <!> which satisfies 
{y^-A^d^}<i>ix) = na-2N^f (4.11) 
where = h“^Vo,V /3 and 

Aa = ^9aN, f = eN~^ dd-iix,xo). (4.12) 

We wish to rewrite this wave-equation in self-adjoint 
form. This ensures that the corresponding Green’s func¬ 
tion will be symmetric in its arguments and hence our 
singular field will exert no force on the charge, but acts 
only to renormalize the moments of the body [27]. We 
can express the electrostatic wave equation as 


N^/hdx°‘\ N dx^) ) 


(4.13) 


which suggests that electrostatic potential satisfies a 
Poisson equation on a conformally related metric with 
a rescaled charge. To see this, we choose our conformal 
metric to satisfy 


Now computing the self-force requires evaluating the 
gradient of the regular field 4 )r = $ — before taking 
the coincident limit x —>■ Xq. Hence only w(xq) = wq 
and its derivative appear in the expression for the self¬ 
force. Notwithstanding the arbitrariness of this function 
which is necessary to encode different boundary condi¬ 
tions in the physical field, one still expects the singular 
field to be insensitive to such boundary conditions and 
hence w(x) for the singular Green’s function ought to be 
locally constructed. We postulate the following axioms 
on w: 


1. w must transform appropriately under 
length rescalings. Explicitly, under a change 
of units of length by a factor I, we have that 


CT ^ l^a, V rW and 
by Eq. (4.17) 


I ^G„ and hence 


w ^ r^(w-vlogl^). (4.20) 

This implies that w assumes the form 

w = V log\Qi\ + Q 2 , (4.21) 


Vhh-'^ = 

N 

which in turn implies 

Kp = 


(4.14) 


(4.15) 


whence Eq. (4.13) may be recast into the Poisson form 

= eHd-2^d-i(®,a:o), (4.16) 

where is the D’Alembertian operator with respect to 
the conformal metric hap and S{x, Xq) = S{x — Xq )/ \/h. 

Restricting; attention to a static charge in a 5D space- 
time, then hap = N~^hap and the Hadamard Green’s 
function corresponding to the operator has the form 


Gh{x,x') = -(^ 


1 f A^/‘^{x, x') 


+ V{x, x') log tr -I- W(a;, a;') 

(4.17) 


where all biscalars appearing here are w ith respect to 
hap- In light of the expansions (4.9 )-(4.10), we have that 

A^/^(a;,a;') = 1 -|- ^Ra>p'^'’“ - ^Ra'p'--y'd-'°‘ a’'* 

+ 0(d^) 

Vo{x,x') = -^R{x') + ^R.y{x')a'’°'' + 0(ct), 

Wo(x,x') = w(x') — -I-0(ct). (4-18) 

The singular field can be obtained by a trivial convolution 
of this c onform al Green’s function Gh against the source 
term in (4.16), 


$s(a;) = -eGH(a;,a:o) 

e / A^/^(a:, ccq) 


-I- V{x, Xq) log ct -I- W(a;, xq 
(4.19) 


where Qi and Q 2 are scalars with dimensions in¬ 
verse length squared [L]~^. 


2. w must be constructed only from local geo¬ 
metrical quantities. This is completely reason¬ 
able since the singular field should only be sensitive 
to the local physics near the particle. Since the sin¬ 
gular field satisfies the Poisson equation (4.16), the 
only geometric scalars in the theory are curvature 
invariants of hap- Since, the only curvature invari¬ 
ant with dimensions is the Ricci scalar on 

hap, the possible terms appearing in Qi and Q 2 are 
of the form F{K) R, where K is a dimensionless ra¬ 
tio of curvature scalars, e.g., RabR°'^/R^, and F is 
an arbitrary dimensionless function of such ratios. 


3. Finally, we assume that w does not depend 
on ratios of curvature scalars. Unlike the pre¬ 
vious two axioms, this is a strong assumption. Ob¬ 
viously we want our prescription to be valid for 
arbitrary static spacetimes and the most straight¬ 
forward way to achieve this is to rule out such 
curvature ratios. The more conservative way to 
achieve this would be to impose analyticity in K 
in some complex neighborhood of the positive real 
axis (see, for example, |36] for a discussion in the 
context of quantum field theory). It turns out 
that this restriction is not strong enough to give 
any predictive power since there is still an infinite 
space of functions that satisfy this analyticity cri¬ 
terion [37], for example, F{K) = Rf^/(1 -|- K^) and 
F{K) = K'^ exp{—K'^) are smooth functions satis¬ 
fying F{K) —>■ 0 as 77 —>■ 0 and as 77 —)■ 00 . This 
axiom implies that 

Qi,2 = 011^2 R^ (4.22) 

where a 1.2 are arbitrary dimensionless constants. 
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Now consideration of the above axioms necessitates that 
w = V logai\R\ + a 2 R, ai > 0. (4-23) 

There is still so me re dundancy in this construction. We 


note from Eq. (4.18) that v{x) = Vq{x,x) = —^R{x), 


and hence we can redefine our constants such that 

w = V log ai|h| + 02 V, oi > 0. (4.24) 

Finally, we can absorb Oi into a redefinition of a^-, for 
example, defining a 2 = a — log a\ results in 


w = V log Ihl + av. 


(4.25) 


The crux of this analysis is that we have reduced the 
arbitrariness contained in the function 'w(x) down to a 
single arbitrary dimensionless constant for an arbitrary 
static 5D spacetime. The constant does not depend on 
the particular spacetime. 

The previous analysis applies almost verbatim for a 
minimally coupled scalar field (or a non-minimally cou¬ 
pled field in a vacuum spacetime) so we omit much of the 
details. The field equation for a static scalar charge in a 
static 5D spacetime can be recast into a Poisson equation 

y'^(p{x) = qNoS{x,Xo), (4.26) 

where now the wave operator is with respect to the metric 
hap = N hap- (4.27) 

The singular field has the form 

^s{x) = ^^’^°'^- + V{x,Xo) logcr 

4 V a 

+W(a:,a;o)y (4.28) 


where all quantities are functions of the metric (4.27). 


The axioms we applied to the electrostatic case in or¬ 
der to reduce the arbitrariness of the coincidence limit of 


Wo{x,x') apply also to the scalar case. Hence, w{ x) has 
precisely the same form as in the electrostatic case (4.25) 
but with 


V 


(4.29) 


where R is now the Ricci scalar of the metric (4.27). 


C. 


Singular Field for a Static Charge in 5D 
Schwarzschild-Tangherlini Spacetime 


We will now derive coordinate expansions for the sin¬ 
gular field of a static electric and scalar charge in the 
5D Schwarzschild-Tangherlini spacetime. For an electric 
charge, the singular field is given by Eq. (4.19) where the 


biscalars therein possess ex pansi ons of the form (4.18) 
and w{x) is given by Eq. (4.25). Moreover, the met¬ 
ric that induces the Poisson form of the electrostatic 
wave equation is hap = f~^^^hap where we recall that 
/ = 1 — r\lr'^ and the Ricci tensor and Ricci scalar on 
this metric are 


Ra H 


R = - 


diag{ - ^^(3/ + 1), ^(3/ - 


2r6y7 


(4.30) 


where ^ab is the metric on the 3-sphere. Combining 
Eqs. (4.18), (4.19), (4.25) and (4.30) with a coordinate 


expansion for a yields a coordinate expansion for the sin¬ 
gular field. If we specialize to the situation where the field 
point X and charge location Xq are separated only in the 
radial direction, then defining Ar = r —rg and expanding 
the expression (4.19) in powers of Ar gives 


„3/2 

4>s(r) = -^^° 


Q 3/2 „l/2 

3e Vq fo' 


2 Ar^ 4r 


Ar 




64r2 


^ i (32 - (29 + 3„)» - 390 log [^]) 


f 

Jo 

(^80- 6(26-b3a)yo + (71-b 15a)yo - 32 /o( 6-5yo)log ) -b O(Ar^) (4.31) 


where /g = /(rg) and we have expressed our expansion 
in terms of an inverse dimensionless radius 


Vo = 


(4.32) 


For the static scalar charge in a 5D Schwarzschild- 
Tang herlin i spacetime, the singular field is given by 
Eq. (4.28) with iVg = -yZ/g. The metric that induces 


the Poisson form of the static scalar wave equation is 
hap = hap for which the Ricci tensor and Ricci 
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scalar are 




diag| 


irl 

2 r 6 / 


2 , 0 , 0 , 0 }, 


(«3) 


As in the electrostatic case, we combine Eqs. (4.181, 
(4.28), (4.25) and ( |4.33 ) with a coordinate expansion for 
a to obtain a coordinate expansion for the singular field, 
which for radial separation yields 


. , 9/0 

5/2 A 

q 2 / 0 ' Ar 


3/2 
9 % 


vl 


Ar 96r2 /g 


(24 - (19 - 3a)yo + 3yo log ) 


64r5 




(16 - 2(10-3a)yo +3(3- 0 ) 2 / 0 + 32 /o( 2 -yo) log ) + 0{Ar'^). (4.34) 

\ LiOrg/gJ/ 


V. ELECTROSTATIC SELF-FORCE 


In order to compute the self-force, we replace the field 


gradient in the formal definition (2.7) with the regular¬ 


ized field gradient before evaluating at the charge’s loca¬ 
tion, 

-P'’ = e ^/fo^r^n{xo) =e^/fo[^r<^ - (5.1) 

Rather than working with the coordinate-dependent self¬ 
force, let us introduce the invariant 

F = ±V 9 abF-F» = /g-'/^’- = e [5,4> - 

(5.2) 


which represents the magnitude of the force measured by 
an observer at the charge’s location. The sign was chosen 
to coincide with the sign of F"^. To compute F requires a 
coordinate expansion for the electrostatic field which we 
derived in closed form in Sec. m Taking the angular 
coincidence limit 7 = 0 in Eq. (2.54) and expanding in 
powers of Ar yields 




3/2 


O 3/2 .1/2 

3e 2/g' fo' 


2 Ar^ Ar„ 


Ar 


e 2/0 

32r2 jV 2 


( - 4 - 282/0 + 352/2 + //^(4 + 62/0) 

log ( ”‘)(;^’ ) H log fjo) - 1 - - ISOilo + 101»o 

8V7o ) y-S/o 


8 v^ ! ' “'"''"“■ri/g. 

a /^(8 -k I 62/0 - 3 O 2 / 0 ) - 62 / 0(6 - 52 /o) log - 32 /o (6 - 62 / 0 ) log (5-3) 


where we have again expressed our series coefficie nts in 
term s of 2/0 = ^h/^o- A dopt ing the expansions (4.31) 
and (5.3) in the definition (5.2) results in the closed-form 


expression for the force invariant (where for convenience, 
we now drop the zero subscript and refer to the particle’s 


coordinate location as r) 


e2 y3/2 


^=6^^1^-8 + 202/ + 6(4-3a-31og 12)2/2 
_15(2 - a - log 12)2/3 + + 8y - 152/2) 

- 1 - 62 / 2(6 — 5y) arctanh ^/f |, 


(5.4) 
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FIG. 2: Plot of the self-force invariant on a static electric charge in a 5D black hole spacetime. We have plotted the force as 
a function of the inverse dimensionless radius y = r\lr^ for a set of a values. We see that the force diverges for all values of a 
as the static charge is placed on the horizon i/ —> 1. However, for one particular value, Ocrit = 2, the leading order divergence 
vanishes and we have a milder singularity on the horizon. We also see that the force is repulsive for a. < Qcrit, while for a > acrit 
the force becomes attractive within some radius that depends on ot. 


where y = r|/r^. We can simplify the form of this ex¬ 
pression by redefining the arbitrary coefficient 

a —l a — log(12) (5.5) 

whence the self-force assumes the form 

2 3/2 

^ = ell ( - 8 + 20y + 6(4 - - 15(2 - a)y^ 

+2y/J{4 + 8y- 15y^) -I- 6y^(6 - 5y) arctanh y/f'j . 

(5.6) 

Evidently, the self-force depends on the arbitrary coef¬ 
ficient a, which was inherited from the arbitrariness of 
Wq{x^Xq) in the Hadamard regularization prescription. 
This coefficient is dimensionless and so there are no un¬ 
resolved length-scales in our expression for the self-force, 
contrary to the results of Ref. m- Moreover, our ex¬ 
pression for the self-force does not depend on any simple 
model for the electric charge and is regular for all r > r^, 
contrary to the calculation of Ref. M- However, our cal¬ 
culation does not inform us whether the first order self¬ 
force in 5D depends on internal structure, since there 
are dimensionless parameters that characterize proper¬ 
ties of a body, and it is conceivable that our variable a 
could depend on such parameters. For example, in 5D 
flat spacetime one such dimensionless parameter is 

^ J d'^xf{x)4>*{x)\x\'^, (5.7) 

which arises from the trace part of the quadrupole cou¬ 
pling of the self-energy density of the body with an ex¬ 
ternal gravitational field. Moreover, such dependence on 


internal structure for spherically symmetric bodies arises 
in the second-order electromagnetic self force in four di¬ 
mensions in flat spacetime |38j . A hint that it may arise 
in the first order self-force in five dimensions is provided 
by the perspective of effective field theory (see Refs. [39l- 
m] for the effective field theory approach to the self-force 
problem). When one writes down all the possible terms 
in an action for a charged point particle, one term that 
is allowed is 

Co J dr (5.8) 

where a is the covariant acceleration, r is proper time 
and e is the charge. (This term will give rise to higher- 
derivative equations of motion but a second order in time 
equation of motion can be obtained by the usual reduc¬ 
tion of order technique). In four dimensions, the param¬ 
eter Co has dimensions of length (in units with c = 1), 
and so the operator is an irrelevant operator: its effect 
becomes negligible for very small bodies. In particular 
Co would be expected to be of order the size of the body. 
However, in five dimensions, cg is dimensionless, and so 
the operator is marginal, and should give a non-vanishing 
contribution to the equations of motion in the limit of 
small bodies. The coefficient cq may be a universal con¬ 
stant, independent of body structure, or it may depend 
on the body’s internal structure. In an upcoming paper 
|42) . we will apply Harte’s formalism pTHSO] to address 
whether the self-force on static bodies in arbitrary dimen¬ 
sions depends on internal structure and further elucidate 
the connection between the ambiguity in the choice of the 
singular field and renormalizations of the body’s multi¬ 
pole moments. 
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The asymptotic form of F at infinity {y —>■ 0) is 

(5-9) 

while near the horizon (y = 1) we have 

/ 2 — a 158 — 63 q; 

128v^rH^^ 4rH (r - r^y/'^ 

+ ^^)+0((^-^h)'/^). (5.10) 


on the ergosphere in the Kerr spacetime [43] since it 
requires infinite acceleration to hold the charge hxed 
there. On the other hand, this “intuition” fails in certain 
cases, for example, the electrostatic self-force invariant 
in Schwarzschild spacetime is F = e^M/r^ [H], which 
is everywhere regular. Presumably, the regularity of the 
force on the horizon is a coincidence of the Schwarzschild 
geometry in four dimensions and there appears to be no 
reason to expect this to be true in general. 

VI. STATIC SCALAR SELF-FORCE 


It is clear that for any choice of a, the force invariant di¬ 
verges as we approach the horizon which is in agreement 
with the near-horizon behavior found in Ref. [14j . We 
note that the leading order pathology can be removed by 
choosing a = 2 but the subleading divergence survives. 
Moreover, the force is everywhere repulsive for a < 2 
but for a > 2, the force remains repulsive far from the 
black hole but is attractive within some critical radius 
that depends on a (see Fig. [^. 

That the self-force diverges as the horizon is ap¬ 
proached seems intuitive since a particle requires infinite 
acceleration to remain static on the horizon. For ex¬ 
ample, the self-force on a static scalar charge diverges 


The scalar self-force is defined by 

F'^ = qfodrmixo) = q fo[drV’ - ( 6 . 1 ) 

The force invariant F is defined by 
F = ±V9abF-F» = 

( 6 . 2 ) 


The radially-separated coordinate expansion for the 
scalar field is obtained by expanding the closed-form rep¬ 
resentation (3.20) in Ar yielding 


3 /2 

2 Ar^ 4rH Ar 32r| /o 


(-4-4yo + 3y2 + /^(4-2yo) 


log (^A^) - log ) + ^ ( - 8 - 4» - 4,2 + „» 

- \//o(8 - 8yo + 6yo) + - 2) log “ 2) log + 0(^^^)- 

^ OV/O ' 


(6.3) 


Substituting Eqs. (4.34) and (6.3) into the definition (6.2) produces a more simple form. 


yields (again dropping the subscript 0 for typographical 
convenience) 


F = 


q2 y3/2 


8 — 12y -I- 6(4 — a — log 4)?/^ 


64r| /3/2 

-(10 - 3a - 31og4)y3 - ^/f{8 - 8y + 

-|-6y^ (2 — y) arctanh ). (6.4) 


Redefining a by 


a a — log 4, 


(6.5) 


F = 


64r3 


3/2 

^ (8 - 12y + 6(4 - a)y^ - (10 - 3a)y^ 


-\/7(8 -8y + 6y^)+ 6y^(2 - y) arctanh a//) . 

( 6 . 6 ) 

As for the electrostatic self-force, the scalar self-force does 
not have any undetermined length-scales and need not 
depend on the internal structure of the charge. 

For large r, the force invariant behaves as 


^^^log(2r/rH)-bO(r ^) 


(6.7) 


Hence, we find the scalar self-force to be repulsive at large 
r which disagrees with the results of Ref. Ell though 
the scaling with r agrees. For the static charge near the 
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FIG. 3: Plot of the self-force invariant on a static scalar charge in a 5D black hole spacetime. We have plotted the force as a 
function of the inverse dimensionless radius y = r^/r^ for a series of a values. As in the electrostatic case, the force diverges 
for all values of a as the static charge is placed on the horizon j/ —> 1, but the divergence is milder for acrit = 10/3. The force 
is repulsive for a < oicrit, while for a > Ocrit, there is a radius which depends on a within which the charge feels an attractive 
self-force. 


horizon, we obtain 

/ 10-3a 78 - 33a \ 

128V2rH^^ ^ (r - / 

+0((r-rH)'/'). (6.8) 

As for the electrostatic case, the force blows up as the 
static charge is placed closer and closer to the horizon 
for every choice of the dimensionless coefficient a. How¬ 
ever, for a = 10/3 the leading order singularity vanishes 
and this value delimits the choices for which the force is 
everywhere repulsive, corresponding to a < 10/3, from 
those for which the force turns over and becomes attrac¬ 
tive as we move towards the horizon, corresponding to 
parameter values a > 10/3 (see Fig. [^. 


VII. CONCLUSIONS AND DISCUSSION 


In this paper, we have computed closed-form expres¬ 
sions for the self-force on static electric and scalar charges 
in the 5D Schwarzschild-Tangherlini spacetime. The cal¬ 
culation was facilitated by two results: First, we de¬ 
rived closed-form representations of the electrostatic and 
scalar static fields. Second, we developed an axiomatic 
approach to constraining the possible forms of the sin¬ 
gular field, reducing the arbitrariness of the Hadamard 
form to a single arbitrary dimensionless coefficient. From 
our closed-form expressions for the self-force, Eqs. (5.6) 


and Eq. (6.6), we immediately deduce that the self-force 


does not depend on any unresolved length scales, such 
as the radius of a sphere centered at the charge’s loca¬ 
tion, nor does our calculation necessitate a dependence 
of the self-force on the internal structure of the charge. 


We note that our regularization scheme is quite robust, 
valid for static charges in arbitrary static 5D spacetimes 
and easily generalizes to higher dimensional static space- 
times. It relies on three axioms, two of which are quite 
innocuous assumptions that the singular Green’s func¬ 
tion scales appropriately under a rescaling of length and 
is constructed entirely from local geometrical quantities. 
The third axiom is a strong constraint that rules out ra¬ 
tios of curvature scalars appearing in the singular Green’s 
function, but it would be surprising if the correct expres¬ 
sion for the self-force obtained by a matched asymptotics 
expansion involved such ratios appearing in some special 
combination that is regular for every static spacetime. 


The strong dependence of the self-force on the parity of 
spacetime dimension is intriguing. It seems intuitive that 
the singular Green’s function should only be sensitive to 
the local physics in the vicinity of the charge and that it 
must be symmetric, and yet there is no obvious analog of 
a Detweiler-Whiting Green’s function in odd dimensions 
satisfying these criteria. Indeed, there are strong indica¬ 
tions that such a Green’s function does not exist. This 
raises the very interesting question as to how to regular¬ 
ize the self-force in odd dimensions. Moreover, it remains 
to be proven whether or not the first order self-force de¬ 
pends on the internal structure of the charge. For static 
configurations, we address these questions in an upcom¬ 
ing paper [32] ■ In the general dynamical case however, it 
may be that resolving some of these issues will require a 
derivation of the self-force using matched asymptotic ex¬ 
pansions, or equivalently, by the method of Gralla, Harte 
and Wald [III- 




















17 


Note 

Upon completion of this work, we were made aware of 
a different derivation by Frolov and Zelnikov of the closed 
form Green’s functions presented here Their 

derivation does not rely on summing the modes as we 
have done but rather relates the static Green’s function 
on the black hole spacetime to the static Green’s func¬ 
tion on the Bertotti-Robinson spacetime by a simultane¬ 
ous rescaling of the induced metric and the lapse. Their 
method is valid for the d-dimensional Schwarzschild- 
Tangherlini or Reissner-Nordstrom spacetimes and in 
general yields a double integral expression for the static 
Green’s function. For d = 5, the integrals can be writ¬ 
ten in terms of Elliptic integral functions and we have 
explicitly checked that our derivation agrees. 
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